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INTRODUCTION 

The  particle  velocity  (u),  vs  time  (t),  behavior  of  detonating 
cast  TNT  has  been  measured  by  the  electromagnetic  velocity  (EMV) 
gage  method.  This  method,  which  was  introduced  into  the  international 
literature  by  Zaitsev,  Pokhil,  and  Shvedov1  in  i960,  is  a  direct 
application  of  Faraday's  law  of  electromagnetic  induction.  The 
particle  velocity  in  a  non-conducting  medium  can  be  obtained  by 
measuring  the  eraf  developed  across  the  base  of  the  EMV  gage  which  is 
moving  with  the  medium  in  a  magnetic  field.  This  magnetic  field  must 
be  oriented  normal  to  the  direction  of  the  gage  base  and  the  direc¬ 
tion  of  motion.  The  generated  emf  across  the  gage  base,  in  volts, 
is  given  by 

V  =  H'U-t.  lCf3  (1) 

where  H  is  the  magnetic  field  in  gauss,  u  is  in  mm/p.sec,  and  l  is 

the  gage  base  length  in  cm.  The  EMV  gage  method  and  the  necessary 

instrumentation  were  described  and  checked  out  in  a  non-conductor 
2  3 

(PMMA).  That  worK  showed  that  the  measured  particle  velocity 
agreed,  within  experimental  error,  with  previous  work  using  a 
different  method  on  PMMA. 

When  used  to  study  the  flow  behind  detonation  waves,  the  EMV 
gage  method,  in  principle,  measures  not  only  the  u, t  behavior  of  an 
explosive  but  can  also  be  used  to  determine  the  Chapman- Jouguet  (CJ) 
particle  velocity  (uCJ)  and  the  reaction  time  (tCJ)  directly  and 

thus  the  reaction  zone  width  and  the  CJ  pressure  (P  ) .  All  of  the 

1  4-11  CJ 

reported  work  *  in  explosives  using  the  EMV  gage  method  is  of 

Russian  origin,  mostly  by  Dremin  or  his  associates.  The  Russian 

results  of  the  CJ  parameters  of  several  explosives  raise  the  question 

of  misinterpretation  of  free  surface  velocity  measurements  by 

4  7  12 

American  researchers.  In  particular  for  cast  TNT,  Dremin  * 

et  al,  have  reported  a  CJ  pressure  which  is  about  13%  lower  than  the 

13 

findings  of  Craig,  who  used  the  free  surface  velocity  method  to 
obtain  the  u~  _  vs  x  curve  for  Dural  on  TNT.  Dremin  has  argued 
that  the  wrong  point  was  selected  in  the  observed  free  surface 
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velocity  vs  distance  curve  as  the  value  to  be  associated  with  the  CJ 
condition.  Craig  argued  that  the  region  between  the  CJ  point  he 
chose  and  the  value  chosen  by  Dremin  is  a  decay  zone  caused  by  the 

14 

divergent  flow  behind  the  detonation  front.  Petrone * s  paper  on 
nitrome thane,  which  is  another  explosive  whose  CJ  parameters  are 
also  in  question,  supports  Dremin1 s  arguments.  Our  results  in  cast 
TNT,  which  are  preliminary,  tend  to  support  Dremin. 

Cast  TNT  was  chosen  instead  of  pressed  TNT  in  this  study  because 
its  reaction  time  should  be  greater  than  the  rise  time  of  our  record¬ 
ing  system,  0.10  usee,  which  is  one-third  of  Dremin’s^  reported 
reaction  time.  This  situation  is  necessary  to  observe  any  change  in 
slope  which  might  be  associated  with  the  CJ  point. ^  The  present 
work  shows  that  a  single  record  cannot  be  used  to  determine  a  CJ 
point  in  cast  TNT;  this  may  be  true  for  other  explosives  as  well. 
Another  approach  to  the  problem,  that  of  using  charges  of  different 
length  and,  thus,  varying  the  steepness  of  the  Taylor  wave,  was  used. 

EXPERIMENTAL 

Setup 

The  explosive  charge  and  booster  configuration  used  in  this  work 
is  shown  in  Figure  1.  A  Baratol-Pentolite  plane  wave  booster, 

PWB  (50*8  or  63,5  rum  diameter),  was  used  in  each  experiment  to  initi¬ 
ate  either  the  cast  TNT  directly  or  an  intermediate  layer  of  a 
different  high  explosive  which  acted  as  a  booster  for  the  TNT.  The 
intermediate  layer,  when  present,  was  either  pressed  Pentolite 
(p0  =  1.56  ±  0.01  gm/cc)  or  pressed  TNT  (pQ  =  1.60  ±  0.01  gm/cc) 
pellets  which  had  been  machined  to  50.8  mm  diameter  and  thickness  L. 
The  particle  size  of  the  TNT  was  ~200  microns.  The  cast  TNT  was 
machined  from  cream  cast  charges  into  pellets  of  50.8  mm  diameter 
and  thickness  x.  (The  TNT  was  melted,  cooled  with  stirring  to  78°C, 
and  then  poured  into  a  15"  long  heated  mold  (75°C).)  The  density  of 
the  cast  TNT  was  1.62  ±  0.01  gm/cc.  The  reported  detonation  velocity, 
D,  of  the  cast  TNT  is  6.85  mm/psec,  of  the  pressed  TNT  is  6.91  mm/psec, 
and  of  the  pressed  Pentolite  is  7.28  mm/psec.  The  intermediate  layer 
was  necessary  because  the  PWB's  we re  not  powerful  enough  to  initiate 
the  cast  TNT  close  to  the  PWB-east  TNT  interface.  Jacobs,  et  al1*^. 
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have  shown  that  cast  TNT  requires  a  run  length  which  is  related  to 
the  initiating  pressure  to  attain  steady  state  detonation;  for  75  kb 
the  distance  is  20mm  and  for  137  kb,  6  mm.  Because  Baratol  would 
induce  pressure  of  about  110  kb,  at  the  interface,  the  run  length  to 
steady  detonation  this  would  produce  in  cast  TNT  would  be  intolerable 
for  our  purposes.  For  this  reason,  an  intermediate  explosive  which 
would  reach  steady  state  detonation  quickly  was  added  (see  Fig.  1). 
Pentolite  and  pressed  TNT  were  chosen  because  their  impedances, pqD, 
are  very  close  to  cast  TNT:  pqD  =  11.35*  11.06  and  11.10  gm-inm/p.sec- 
cc  for  Pentolite,  pressed  TNT,  and  cast  TNT  respectively. 

The  PWB  was  initiated  by  a  primacord  lead  (120  grain  EDX/foot) 

30  cm  long  which  in  turn  was  initiated  by  an  exploding  bridgewire 
detonator.  The  primacord  isolates  the  charge  from  the  detonator  to 
prevent  possible  stray  electrical  signals  from  the  firing  unit  being 
picked  up  by  the  gage. 


The  usual  EMV  gage  consists  of  a  rectangular  loop  of  aluminum 
foil  0.13  mm  (5  mils)  thick,  5  mm  wide,  and  10  mm  long.  It  is 
mounted  in  a  cast  TNT  back-up  assembly  whose  thickness,  F,  is  25.5  mm 
(see  Fig.  1).  The  gage  is  mounted  in  the  cast  TNT  by  shaping  the 
foil  around  piece  B.  A  thin  layer  of  silicone  grease  is  placed 
between  the  pieces  except  near  the  gage.  Then  pieces  B,  C,  and  D 
are  placed  together  and  cemented  using  Duco  on  the  exterior  under  a 
slight  pressure.  The  gage  circuit  is  completed  by  connecting  the 
foil  leads  (30-55  mm  length)  to  a  RG  58  C/U  coaxial  line  (50  ohms 
nominal  impedance)  with  a  resistor,  R  (50  ohms),  in  series  with  the 
foil. 


A  0.05  mm  thick,  6.4  cm  square  aluminum  baffle  (see  Fig.  l)  was 
located  in  the  plane  of  the  HE-cast  TNT  interface  and  was  grounded . 
The  significance  of  this  baffle  will  be  discussed  in  the  following 
section. 


The  magnets 
discussed  in  Ref 
Rawson-Lush  Type 


and  power  supplies  used  in  these  experiments  are 
(2)  and  (3).  The  magnetic  field  was  measured  with  a 
824  rotating  coil  gaussmeter  (+  1$  accuracy)  at  the 
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mean  position  of  the  gage.  Then,  with  no  field  present,  the  charge 
assembly  was  placed  in  position,  the  current  and  therefore  the  mag¬ 
netic  field  returned  to  the  same  value  as  before  and  the  shot  fired. 

Instrumentation 

The  EMV  gage  was  connected  to  15  feet  of  RG  58C/U  (50  ohm 
nominal  impedance)  coaxial  cable  with  a  50.8  +  0.25  ohm  resistor  in 
series  with  the  gage  to  minimize  any  reflections.  The  cable  was 
connected  to  the  center  lead  of  a  50  ohm  power  divider  (Tektronix 
Model  017-0082-00).  The  other  two  leads  of  the  power  divider  were 
connected  to  two  fast  response  oscilloscopes,  a  Tektronix  Model  454 
(2.4  ns  rise  time)  and  a  Hewlett-Packard  Model  175  (7*0  ns  rise  time), 
through  50  ohms  terminators  (Tektronix  Model  011-049-01).  Figure  2 
is  the  circuit  diagram  of  this  arrangement.  The  arrangement  results 
in  the  prcper  termination  of  all  cables;  that  is.  the  coaxial  cables 
see  50  ohms  at  both  ends.  This  arrangement,  however,  results  in  a 
considerable  drop  in  amplitude  of  the  signal  generated  by  the  gage. 
The  50  ohm  resistor  near  the  gage  reduces  the  signal  by  50#  while  the 
power  divider  further  splits  this  signal  by  54.88#  with  the  net 
result  that  each  oscilloscope  sees  only  27.44#  of  the  original 
voltage . 

Voltage  calibration  marks  were  placed  on  all  oscilloscope  camera 
records.  This  voltage  was  obtained  from  a  power  supply,  the  output 
of  which  was  measured  with  a  Fairchild  Model  7050  Digital  Multimeter 
(+0.1#  accuracy).  The  calibration  voltage  was  placed  on  the 
oscilloscope  records  with  the  power  divider  in  place;  this  procedure 
eliminates  any  corrections  due  to  imperfect  splitting  of  the  signal 
by  the  power  divider.  Time  calibration  was  obtained  from  a  crystal- 
controlled  time  mark  generator  (Tektronix  Model  180A)  and  was  placed 
on  separate  records  because  of  the  limited  viewing  area  of  the 
oscilloscope  screen.  The  oscilloscopes  were  triggered  internally  by 
the  output  voltage  from  the  gage. 

The  aluminum  baffle  shown  in  Figure  1  was  found  useful  in 
reducing  spurious  signals  which  appeared  in  the  first  experiments 
with  cast  TNT;  this  phenomenon  also  appeared  in  the  work  with  PMMA2* ^ 
where  this  baffle  was  also  employed  for  the  same  purpose.  (Figure  5 
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POWER  DIVIDER 


FIG.  2  CIRCUIT  DIAGRAM 
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HOT  REPRODUCIBLE 


k 


1CH  SHOT  130 


T 

0.6  VOLTS 

_L 


GAGE  LOCATED  5C  .8  MM  FROM  PENTOLITE  BOOSTER 

MAGNETIC  FIELD  750  GAUSS 
SWEEP  SPEED  0.05  iiSEC/CM 


FIG.  3  EXAMPLE  OF  A  RECORD  OBTAINED  WHEN  NO  ALUMINUM  BAFFLE  WAS  USED. 
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shov;s  an  example  of  the  signal  obtained  without  a  baffle.)  The  baffle 

was  grounded  to  a  circuit  which  was  separate  from  the  oscilloscope 

circuit.  This  grounded  baffle  greatly  reduced  and,  in  most  cases, 

completely  eliminated  the  spurious  signals.  A  similar  type  of  noise 

T  6 

was  observed  by  A.  A.  Brisn,  et  al,"  when  the  electromagnetic  method 
was  used  to  study  conductivity  in  detonating  explosives. 

» 

Data  Reduction 

The  oscilloscope  records  were  read  on  a  Universal  Telereader  with  v 

the  origin  at  the  point  where  the  trace  leaves  the  baseline.  The 
digital  data  obtained  were  then  converted  to  u  and  t  by  means  of  a 
computer  program  using  eq.  (1)  and  the  calibration  of  the  oscillo¬ 
scope  . 

FACTORS  WHICH  AFFECT  RECORDS  AM?  THEIR  INTERPRETATION 

There  are  several  factors  which  must  be  considered  when  inter¬ 
preting  the  records.  These  can  be  divided  into  two  main  groups: 

(a)  those  which  contribute  to  the  rise  time,  and  (fc)  those  which  are 
connected  with  the  conduction  of  the  plasma. 

System  Rise  Time 

There  are  at  least  five  factors  which  contribute  to  the  overall 
rise  time  of  the  recording  system;  they  are: 

1.  wave  front  curvature, 

2.  shock  impedance  mismatch  between  the  gage  and  the  cast  TNT, 

3.  'wave  front  tilt, 

4.  Inherent  rise  time  of  the  oscilloscopes, 

5.  electrical  impedance  mismatch  between  the  coaxial  cable  and 
the  oscilloscope  or  the  EMV  gage.  The  first  three  are  connected  with 
the  detonation  itself,  while  the  last  two  are  electrical.  These 
effects  have  been  discussed  in  Ref.  (2)  and  (3).  The  planarity  of  the 

wave  front  after  it  had  traveled  from  a  50.8  mm  diameter  PWB  through  v 

25.4  mm  of  pressed  TNT  and  25.4  mm  of  cast  TNT  was  examined  using  a 
streak  camera  (writing  speed,  4  mm/ixsec);  the  front  was  plane  to 
within  reading  error  over  the  middle  20  mm  of  the  charge  which  allows 
factor  1  to  be  neglected,  see  Figure  4.  However,  a  tilt  of  0.2  mm,  or 
0.02  radians  was  observed  in  the  region  where  the  gage  would  have 
been  placed. 
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This  would  result  in  a  rise  time  of  29  ns  for  factor  3>  T^.  The 
total  rise  time,  T  ,  can  be  calculated  using  the  equation 

To  -Jf  Ti2"  <2> 

p  ^ 

where  the  Ti ’s  are  the  individual  contributions  listed  above.  * 

If  we  let  T^  =  0  and  29  ns,  Tq  is  87  and  92  ns  for  T^  =  87  ns  (two 
double  reflections  in  the  pickup  foil  using  an  average  shock  velocity 
of  6  mm/p.sec).  If  three  double  reflections  are  used,  T  oecomes  130 
and  134  ns  respectively.  Table  1  lists  the  time  observed 


for  the  trace  to  reach  the  maximum  point  in  the  curve,  Tmax,  for 
experiments  in  which  a  5  mil  fell  was  used.  Tq  was  estimated  by 
taking  90#  of  Tmax.  All  but  three  of  the  experiments  have  an 
estimated  rise  time  below  134  ns.  The  longer  rise  times  are  attri¬ 
buted  to  a  tilt  greater  than  that  shown  in  Figure  4.  Thus,  wave 
front  tilt  can  have  an  appreciable  effect  on  T  ,  and  it  can  be  larger 
than  the  specific  example  given. 


The  effect  of  rise  time  on  the  true  signal  was  investigated 

p  -3 

theoretically  *  .  A  triangular  pulse  whose  initial  fall  off  was 
rapid  and  then  abruptly  became  less  steep  (i.e.,  double  sloped)  was 
used  as  the  input  to  an  RC  circuit.  A  double  sloped  pulse  was  chosen 
because  it  is  an  approximation  to  the  von  Neumann  spike  followed  by 
a  Taylor  wave.  This  calculation  showed  that  the  observed  trace  will 
overshoot  the  applied  pulse  at  some  time  after  the  driving  pulse 
has  started  to  decay  and  then  eventually  closely  parallel  the  applied 


pulse.  The  problem  was  made  dimensionless  by  introducing  the  vari¬ 
ables  y  =  t/t^,  0  =  T/t^,  and  Z1  =  ei/eQ1  where  ^  is  the  time  to 
the  break  in  the  input  pulse,  t  =  RC  is  the  time  constant  of  the 


circuit,  and  eQ1  is  the  amplitude  of  the  pulse  at  the  break.  It  is 

shawn in  Ref.  2  that  T  ~  2.2t.  The  results  of  this  calculation  are 

o 

shown  in  Figure  5*  The  calculation  also  showed  that  a  recording 
system  with  a  rise  time  smaller  than  0.44  of  the  time  to  the  break 


is  needed  to  locate  the  break.  If  the  reaction  time  in  cast  TNT  is 

4 

O.30  iisec  as  Dremin  claims  then  our  recording  system  will  generally 
be  able  to  reproduce  this  break. 
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Conductivity 

The  effect  of  conduction  is  very  pronounced  in  detonating  cast 
TOT?  and  is  observable  in  two  ways :  (a)  different  output  voltage-time 
curves  obtained  when  the  base  is  left  off  the  gage  and(b)  different 
u,  t  curves  obtained  for  different  A1  gage  thicknesses.  Three  experi¬ 
ments  were  attempted  in  which  the  base  of  the  gage  was  removed. 

Figure  6  is  a  comparison  of  the  results  obtained  for  a  baseless  and 
a  regular  gage  in  cast  TNT.  The  pseudo  particle  velocity  obtained 
from  the  baseless  gage  lies  0.5  mm/p.sec  (~33^)  below  the  particle 
velocity  obtained  from  the  regular  gage.  In  an  attempt  to  gain  more 
information,  the  cable  of  the  system  was  changed  from  50  ohms  (with 
a  50  ohm  resistor  in  series  with  the  gage)  to  93  ohms  (with  a  93 
ohm  resistor  in  series  with  the  gage).  The  results  of  this  experiment 
are  ~0.9  mm/jisec  (~60$)  below  the  regular  gage  valve;  it  was 
anticipated  that  the  curve  would  lie  between  the  regular  gage  and 
the  50  ohm  baseless  results.  The  full  meaning  of  these  results  are 
not  understood  at  present  and  further  work  is  planned  in  this  area. 

In  surprising  contrast  to  the  present  results  from  the  detonation  of 

east  TNT,  the  absence  of  the  gage  base  in  detonating  pressed  tetryl 

17 

did  not  greatly  affect  the  shape  of  the  output-time  curves. 

Dremin^  has  used  various  thicknesses  of  gages  made  of  several 
different  metals.  He  reported  that  gage  thickness  did  not  affect 
the  output,  provided  the  thickness  was  less  than  6  mils.  In 
attempting  a  spot  check  of  this  statement  we  found  that  in  cast  TOT, 
gage  thickness  does  affect  the  results.  Several  experiments  were 
tried  with  A1  and  copper  foils  (resistivity,  2.85  x  10  ^  and  1.73  x 
10~^  ohm-cm  respectively)  of  different  thickness.  Figure  7  shows 
the  results  for  1  and  5  mil  A1  gages  where  all  other  factors  are 
the  same.  As  is  evident  from  the  figure,  the  thinner  and,  therefore, 
higher  resistance  gage  gives  much  lower  values  of  u,  and  u^j  much 
lower  than  any  generally  accepted  values  for  TOT.  Figure  8  shows 
the  results  of  a  5  niil  copper  gage  and  a  5  mil  A1  gage  in  otherwise 
identical  experiments.  (There  is  some  spurious  signal  on  the  copper 
curve;  this  shot,  however,  was  fired  without  the  grounded  A1  baffle.) 
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FIG.  6  RESULTS  OF  REGULAR  AND  BASELESS  GAGE 
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FIG.  7  RESULTS  USING  DIFFERENT  THICKNESS  AL  FOILS 
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Figure  8  also  shows  the  results  for  a  5  mil  copper  gage.  The  curves  of 
the  3  mil  Cu  and  5  mil  A1  in  Figure  8  are  coincident  to  3^.  On  the 
basis  of  these  results,  the  5  mil  A1  gage  was  selected  for  the 
present  work. 

APPLICATION  OF  THE  EMV  METHOD  AND  REPRODUCIBILITY 

A  total  of  24  experiments  have  been  carried  out  in  cast  TNT. 

Of  these  only  one  experiment  resulted  in  no  record  being  obtained. 

Four  of  the  remaining  23  had  a  spurious  signal  superimposed  on  the 
trace  but  still  produced  some  useful  information.  These  experiments 
are  listed  in  Table  2. 

Reproduciblity 

There  are  two  types  of  reproduciblity  involved  in  this  work: 

(a)  that  of  the  u,  t  curves  obtained  from  the  two  oscilloscopes  for 
the  same  experiment;  and  (b)  that  of  the  u,  t  curves  obtained  for 
replicate  experiments.  Figure  9  shows  the  u,  t  curves  obtained  from 
the  two  oscilloscopes  for  shot  138;  the  curves  are  coincident  to 
within  reading  error  (+1$)  for  the  time  covered  by  the  oscilloscope 
using  the  faster  sweep.  The  oscilloscope  with  the  faster  sweep  was 
able  to  record  the  rise  of  the  signal  whereas  the  other  oscilloscope 
did  not  because  of  the  slower  sweep  speed.  The  assumption  is  made 
that  the  records  would  be  coincident  for  the  entire  range;  the  fact 
that  the  u,t  curves  also  agree  for  replicate  shots  (see  below)  adds 
weight  to  this  assumption.  Figure  10  shows  the  u,t  curves  obtained 
from  replicate  shots  (x  +  l  =  38. 1  mm,  shots  138  and  144);  the  two 
curves,  obtained  from  the  HP  175  oscilloscope,  agree  to  within  +1.5% 
for  the  region  100  <  t  <  500  us  and  to  within  +4%  over  the  remainder 
of  the  curves.  The  difference  in  the  first  100  ns  is  due  to  differ¬ 
ences  in  the  rise  time  from  shot  to  shot.  Thus,  reproduciblity  is 
quite  good. 

Evaluation  of  Single  Record  Results 

As  mentioned  in  the  Instrumentation  Section,  the  oscilloscopes 
were  triggered  internally;  the  signal  from  the  EMV  gage  Is  transmitted 
directly  to  the  trigger  circuits  of  the  oscilloscopes  but  is  delayed 
in  transmission  to  the  oscilloscopes  deflection  plates.  Figure  11  is 
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FIG.  8  COMPARISON  OF  COPPER  AND  ALUMINUM  GAGES 
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a  set  of  typical  records  obtained  from  a  single  experiment;  these 
records  correspond  to  the  u,t  curves  plotted  in  Figure  9-  As  can  be 
seen  from  Figures  3,  10,  and  11,  no  sharp  break  or  sharp  change  of 
slope  is  evident  on  the  records.  This  is  true  for  all  but  two  records 
for  cast  TKT.  Figure  12  shows  an  example  of  a  record  with  a  break 
(shot  174).  The  undulations  in  the  trace  over  the  last  half  of  the 
record  are  probably  die  to  motion  of  the  leads  (which  are  connected 
to  the  cable)  because  the  undulations  occur  after  the  detonation  has 

reached  the  end  of  the  charge.  The  infrequent  appearance  of  a  sharp 
break  leads  to  the  conclusion  that  the  CJ  condition  cannot  be  obtained 
from  a  single  experiment.  The  implication  of  similar  work  in  the 

h. 

Russian  literature  was  that  the  CJ  point  could  be  determined  from  a 
single  experiment  if  the  rise  time  of  the  system  was  shorter  than  the 
reaction  time.  Dremin ! 3  records  of  cast  TNT  (Figs.  2a  and  2b  of  Ref.  4) 
do  indeed  show  a  sharp  break.  'When  these  records  are  plotted  to  the 
time  scale  used  in  this  work,  the  sharp  break  disappears;  Dremin’ s 
curves  also  agree  in  shape  with  ours  (see  Fig.  13).  Because  the  break 
in  the  curve  can  be  enhanced  or  reduced  by  using  different  time 
scales  for  the  u,t  plot,  other  approaches  are  desirable  to  insure  that 
the  break  is  meaningful.  One  such  approach  is  described  in  the  next 
section. 

Comparison  of  Observed  and  Predicted  u. t  Curves 

The  discussion  in  Evaluation  of  Single  Record  Results  section 

pointed  out  the  fact  that  distinct  breaks  in  a  single  record  were 

infrequently  observed  and,  therefore,  u  and  the  reaction  time  could 

l/J 

not  be  determined  with  the  EMV  gage  method  using  a  single  record  with 
any  degree  of  certain! ty.  Another  approach,  which,  in  principle,  can 
be  used  to  determine  uCj,  is  to  place  EMV  gages  at  various  distances 
from  the  boosterlng  system  and  compare  the  u, t  curves  obtained  for 
different  lengths  of  run.  This  approach  has  been  used  by  the  Russians^ 
in  a  situation  where  they  could  not  observe  a  break  which  could  be 
associated  with  the  end  of  the  reaction  zone.  The  u, t  curves  should 
coincide  in  the  reaction  zone  region  regardless  of  the  length  of  the 
charge,  because  the  reaction  zone  is  assumed  to  be  independent  of  the 
charge  length  if  the  detonation  has  reached  steady  stace.  The  falloff 
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FIG.  11  TYPICAL  RECORDS 
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FIG.  12  RECORD  SHOWING  DISTINCT  BREAK 
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in  the  u,t  curves  from  the  CJ  value  is  called  the  Taylor  wave^in 
which  flow  variables  depend  on  the  run  length,  x  -f  L,  of  the  charge. 
That  is,  as  x  +  L  is  increased  the  slope  of  the  u, t  curve  will 
decrease  if  the  flov;  is  one-dimensional.  Thus,  by  plotting  the  u,t 
curves  with  the  origin  of  the  plot  coinciding  with  the  zero  time  of 
the  individual  records,  one  should  observe  a  region  where  the  curves 
coincide  and  then  a  region  in  which  they  diverge. 

Qualitatively,  the  u, t  curves  behave  as  mentioned  above. 

Figure  14  is  a  plot  of  the  u,t  curves  obtained  for  charges  having 
(x  +  L)/d  <  1.  For  such  charges,  the  flow  should  still  be  one¬ 
dimensional  and  the  slope  of  the  u, t  curves  should  decrease  with 
increasing  x  +  L.  This  is  the  situation  shown  in  Figure  14.  The 
curves  coincide  for  t  <  300  ns  and  diverge  from  t  ~  300  ns.  Figure  15 
is  a  plot  of  the  u, t  curves  obtained  for  (x  +  L)/d  >1.  In  this 
region,  the  flow  probably  is  no  longer  one-dimensional  because 
lateral  rarefactions  have  had  time  to  reach  the  center  of  the  charge. 
The  one -dimensional  behavior  of  the  u, t  curves  in  the  region  follow¬ 
ing  the  reaction  zone  will  no  longer  be  present.  Figure  15  shows  that 
the  region  of  coincidence  for  t  <  300  ns  is  still  the  same,  but  the 
slopes  of  the  u, t  curves  for  t  >  300  ns  are  not  as  reproducible  and 
do  not  exhibit  a  trend  with  increase  in  (x  +  L).  Thus,  by 
qualitatively  comparing  the  u,  t  curves  for  charges  of  different 
length,  one  arrives  at  a  CJ  velocity  of  ~1.60  mm/usec  and  a  reaction 

time  of  ~  300  ns  which  are  in  excellent  agreement  with  Dremin's 
4 

results  . 


The  experimental  curves  discussed  above  can  be  compared  w ith 
computed  curves  for  plane  one-dimensional  flow  in  a  polytropic  gas. 
The  flow  is  do  scribed  by 


u  = 


2(D-Ugj) 

Y-i 


O) 


where 

u  «=  particle  velocity 

UCJ=  particle  velocity  at  the  CJ  point 

Y  *=  adiabatic  exponent  of  the  gas 
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tR  =  (x  +  L)/D,  the  time  required  for  the  detonation  to  move 
from  the  PWB  to  the  gage 
t  =  reaction  time,  tCJ 
D  *=  detonation  velocity,  6.85  mm/nsec 
\i  =  {y  -  1)  /  (7  +  1) 
x  =  length  of  the  cast  TM?  specimen 
L  =  length  of  the  pressed  TMF  booster 
Equation  3  is  the  relation  between  u  and  t  in  the  Taylor  wave  which 
follows  the  reaction  zone  in  the  simple  model  of  the  detonation  process. 
Some  comments  on  the  derivation  of  the  relation  are  given  in  Appendix 
A.  We  force  the  computed  curve  to  fit  the  experimental  curve  at  the 
break  by  using  the  value  of  u  at  the  break  as  the  value  of  u__  We 

vu  • 

also  use  the  reaction  time,  t,  and  the  total  length  of  the  charge, 
x  +  L,  to  evaluate  This  latter  assumption  cannot  be  correct;  the 
detonation  starts  in  the  PWB.  Because  the  flow  in  the  PWB  is  two- 
dimensional,  there  really  is  no  point  (or  plane)  in  the  PWB  at  which 
both  the  detonation  and  the  following  rarefaction  wave  originate. 

Hence,  we  can  use  this  model  for  only  limited  guidance. 


Figure  16  is  a  plot  of  the  results  from  Shot  161  for  which 

x  +  L  =  38.I  mm.  Computed  Curve  1  is  based  onuAT  =  1.55  nsec  and 

- 

t  =  0.60  nsec.  It  coincides  with  the  gage  record  in  the  time  interval 

0.60  <  t  <  1.2  p.sec  and  parallels  it  for  1.3  <  t  <  1.7  nsec.  Curve  2 

is  based  onu^j  =  1.6  mm/nsec  and  r  =  0.30  nsec.  This  curve  lies 

below  the  experimental  curve  and  parallels  it  for  t  >  0.6  nsec.  If 

we  relied  on  these  comparisons  to  pick  the  values  of  u  and  t  for 

CJ 

Shot  161,  we  would  take  the  values  used  for  Curve  1,  1.55  mm/nsec  and 
0.60  nsec.  This  would  be  risky  because  of  the  uncertainty  in  the 
value  we  have  used  for  the  run  length,  x  +  L,  as  mentioned  above. 


■3 

It  has  been  suggested^  that  the  longer  reaction  time  is  more 
appropriate.  This  is  justified  by  assuming  that  the  TNT  is  not 
totally  reacted  after  0.3  nsec.  Perhaps  it  is  99$  reacted  after  this 
time,  the  other  1 %  reacting  during  the  next  0.3  nsec,  giving  the 
level  portion  of  the  experimental  curve. 
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COMPARISON  OF  PRESENT  RESULTS  WITH  PREVIOUS  WORK 

4 

The  EMV  gage  method  was  used  by  Dremin  and  Shvedov  to  obtain  CJ 
values  for  cast  TNT.  Their  two  records  (Fig*  2a  and  2b,  Ref'  4)  are 
plotted  in  Figure  13  along  with  shot  l6l.  Their  curves  have  been 
displaced  0.1  and  0.2  mra/psec  for  comparison  of  detail;  the  gaps  in 
Dremin* s  curves  are  due  to  the  superposition  of  time  marks  which 
obscure  the  record.  Curve  2  (premin's  Fig* 2a)  and  curve  3  (shot  161)) 
both  have  a  level  region  for  300  <  t  <  600  ns.  Both  of  Dremin* s 
curves  roughly  parallel  shot  161.  The  fact  that  the  curves  are  par¬ 
allel  even  though  run  lengths  are  different  (Dremin  *s  records  are 

for  (x  +  L)  >  2.5  d)  is  accidental;  the  fact  that  level  regions  are 

“  12 
present  and  ignored  by  Dremin  is  significant.  Craig  has 

investigated  pressed  TNT  at  a  density  very  close  to  that  of  cast  TNT. 

He  used  the  free  surface  method  with  Dural  as  the  inert  and  observed 

two  breaks  in  his  u~  vs  x  curve .  He  chose  the  first  break  as 

X  •  s  • 

being  associated  with  the  CJ  point  while  the  region  between  the  two 

breaks  was  described  as  being  a  "decay"  zone.  Veretennikov,  Dremin, 

1  o 

et  al,  have  shown  that  the  double  break  observed  by  Craig  can  be 
eliminated  by  using  a  spall  plane;  this  results  in  the  elimination  of 
the  first  break.  Table  3  lists  our  results  along  with  Dremin* s  and 
Craig's.  As  can  be  seen  from  this  table,  we  are  in  good  agreement 
with  Dremin  and  with  Craig’s  value  of  pressure  at  the  second  break 
point.  Craig’s  first  break  point  results  in  a  CJ  pressure  which  is 
15$  higher  than  Dremin *s  and  ours.  The  accumulation  of  evidence 
(including  the  computation  by  Petrone^)  makes  it  probable  that  Craig’s 
interpretation  of  his  data  was  incorrect. 

CONCLUSIONS  FOR  CAST  TNT 

Reproduciblity  of  records  for  cast  TNT  from  charge  to  charge  is 
within  +3$  for  times  of  a  few  multiples  of  the  reaction  time,  tCJ. 

An  aluminum  gage  having  a  thickness  of  5  mils  (0.13  mm)  seems 
satisfactory  for  work  in  cast  TNT.  Thicker  gages  might  extinguish  the 
detonation  and  would  lengthen  the  rise  time  as  explained  in  System 
Rise  Time  section.  Gages  1  mil  (0.03  mm)  thick  give  results  which 
are  too  low. 
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The  CJ  parameters  of  cast  TNT 

UCJ  =  1-60 
and  tCj  =  0.30 

These  values  are  in  agreement  with 
break  in  Craig's  data1^. 


determined  in  this  study  are: 

ram/iisec 

jisec 

h 

Russian  work  and  with  the  second 


Phenomena  not  yet  completely  understood  are  the  records  obtained 
with  the  baseless  probes  and  the  effect  of  gage  thickness  on  the 
output.  These  have  been  tentatively  attributed  to  the  effect  of 
conductivity  of  the  detonation  products  in  the  first  case  and  of 
the  relative  conductivity  of  the  products  as  compared  to  the  gage  in 
the  second  case. 


GENERAL  CONCLUSIONS 

A  single  record  cannot  be  relied  on  to  give  uCJ  and  tCj. 
Comparison  of  records  for  which  the  detonation  front  has  traveled 
different  distances  removes  most  of  the  ambiguity  in  determining 

UCJ  311(1  tCJ* 

When  applying  this  method  to  other  explosives,  exploratory  work 
should  be  carried  out  at  two  or  more  gage  thicknesses. 
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Table  1 

Estimated  Rise  Times 


Shot 

ji 

V 

Observed 

T  „  (ns) 
max'  1 

EST. 

T  (ns) 
o' 

125 

108 

97 

126 

72 

65 

129 

75 

68 

130 

55 

50 

134 

171 

154* 

138 

45 

41 

144 

50 

45 

143 

69 

62 

147 

52 

47 

148 

70 

63 

153 

218 

175 

155 

84 

76 

161 

133 

120 

163 

168 

151 

164 

98 

88 

166 

146 

131 

174 

88 

79 

175 

108 

97 

*Copper  5  mil  foil 3  all  others  are  A1  5  mil  foil. 
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Table  2 


Individual  Experiments 


X  -rL 

X 

L 

Shot  Inter- 

Foil  Magnetic 

PWB 

(mm) 

^mm) 

(mm) 

No.  mediate 

HE* 

Material  Field 

(gauss) 

Dia. 

(mm) 

Baseless  Gage 


31.8 

25.4 

6.4 

168 

— 

A1 

760 

63.5 

31.8 

2^.4 

6.4 

170 

PT 

A1 

895 

63.5 

76.2 

76.2 

0.0 

128 

FT 

A1 

845 

50.8 

1.0  mil  Gage 

50.8 

5C.S 

0.0 

123 

A1 

835 

50.8 

76.2 

76.2 

0.0 

124 

- 

A1 

788 

50.8 

3.0  mil  Gage 

63.5 

50.8 

12.7 

131 

P 

Cu 

790 

50.8 

5  mil  Gage 

12.7 

6.4 

6.4 

155 

PT 

A1 

890 

63.5 

19.1 

12.7 

6.4 

164 

PT 

Al 

435 

63.5 

25.4 

12.7 

12.7 

153 

PT 

A1 

870 

63.5 

25.4 

12.7 

12.7 

163 

PT 

Al 

642 

63.5 

38.1 

25.4 

12.7 

138 

p 

Al 

255 

50.8 

38.1 

25.4 

12.7 

144 

P 

Al 

258 

50.8 

38.1 

25.4 

12.7 

l6l 

PT 

Al 

840 

63.5 

38.1 

25.4 

12.7 

174 

P 

Al 

868 

63.5 

38.1 

25.4 

12.7 

166 

PT 

Al 

847 

63.5 

50.8 

50.8 

0.0 

126 

- 

Al 

775 

50.8 

63.5 

50.8 

12.7 

147 

P 

Al 

800 

50.8 

63.5 

50.8 

12.7 

130 

P 

Al 

750 

50.8 

76.2 

50.8 

25.4 

148 

PT 

Al 

750 

50.8 

76.2 

76.2 

0.0 

125 

— 

Al 

800 

50.8 

76.2 

76.2 

0.0 

145 

— 

Al 

256 

50.8 

76.2 

76.2 

0.0 

175 

— 

Al 

839 

63.5 

88.9 

76.2 

12.7 

13^ 

P 

Cu 

792 

50.8 

*P  =  Pentolite 
PT  =  Pressed  TNT 


Comparison  of  Results  for  TNT 
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APPENDIX  A 

Derivation  of  One-Dimensional  Plane  Flow  Equation 


Equation  (3)  can  be  derived  by  starting  with  equation  (46. 06) 
from  Courant  and  Friedrichs*.  This  equation  is 


x  =  -(ii-2-l)e0t  +  ^CotoQ-)1-^ 

(Al) 

which  relates  distance  end  time  along  the  path  of  a  particle 
in  an  adiabatic  expansion  process  and  where 

x  =  distance, 
t  =  time, 

t0  =  time  for  detonation  to 
reach  particle, 

c0  =  sound  speed  ahead  of  wave. 

We  also  have  the 

following  relations: 

n 

(A2) 

or 

..-2  _  V+l 

11  "  y=t 

(A3) 

and 

-2  2 

H  -1  =  — T 

y-1  . 

(A4) 

We  now  substitute  equation  (A3)  and  (A4)  into  equation  (Al)  to 
obtain: 

2 

X  =  “  y=I  Cot  +  Cot°(l7)  .  (A5) 

To  obtain  the  equation  for  u(t)  we  differentiate  equation  (A5)  with 
respect  to  t.  This  yields 


*R.  Courant  and  K.  6.  Freidrichs,  "Supersonic  Flow  and  Shock 
Waves" ,  Interscience  Publishers,  Inc.,  New  York,  1948,  p.  104. 
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In  order  for  this  expression  to  represent  the  flow  in  the  Taylor  wave, 
boundary  conditions  must  be  imposed  at  t=t0=tD  +  t  ,  at  which  time 
u  =  uCj.  tR  is  the  time  required  for  the  detonation  to  reach  the 
particle  while  t__  is  the  reaction  time.  The  final  equation  is 


(A7) 


In  the  computations  u.T  and  D  are  the  input  parameters  while  c  and 

C /o  ^ 

y  are  calculated  from 


and 


Equation  A7  is  evaluated  for  t  £  tR  +  t. 


